Cooperativity Beyond Caging: Generalized Mode Coupling Theory 
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The validity of mode coupling theory (MCT) is restricted by an uncontrolled factorization approx- 
imation of density correlations. The factorization can be delayed and ultimately avoided, however, 
by explicitly including higher order correlations. We explore this approach within a microscopically 
motivated schematic model. Analytic tractability allows us to discuss in great detail the impact of 
factorization at arbitrary order, including the limit of avoided factorization. Our results indicate 
a coherent picture for the capabilities as well as limitations of MCT. Moreover, including higher 
order correlations systematically defers the transition and ultimately restores ergodicity. Power-law 
divergence of the relaxation time is then replaced by continuous but exponential growth. 

PACS numbers: 64.70.Pf, 61.20.Lc, 82.70.Dd, 61.20.-p 



The dynamic scaling in dense supercooled liquids and 
colloidal suspensions is a subject surrounded by contro- 
versy. Mode coupling theory (MCT) has shaped our un- 
derstanding through detailed and remarkably successful 
predictions [l| . A salient result of MCT is the existence 
of an ideal glass transition. The latter, however, is pre- 
dicted to occur at substantially higher temperatures or 
lower densities than that observed in the laboratory 
Furthermore it has been suggested that there are "acti- 
vated processes" not accounted for by MCT which re- 
store ergodicity and thus round off the ideal transition 
0, 0] . In order to remedy such shortcomings one should 
scrutinize the approximations made within MCT. 

A fundamental quantity in MCT is the two-point cor- 
relator Fk(t) = A/" -1 (p-k{0)pk(t)) with pk the Fourier 
transform of spatial particle density fluctuations. For J\f 
particles interacting through a pair potential one finds, 
using standard projection operator techniques 



F k {t) + PkF k {t) + ( 7fe * F k )(t) = 0. 



(1) 



Here (/ * g)(t) — L drf(r)g(t — r) and dots indi- 
cate time-derivatives. Further pk — ksTk 2 /mSk where 
ks is the Bolzmann constant, T the temperature, m 
the particle mass and Sk — Fk(Q) the static struc- 
ture factor. The memory function 7fc(£) is related to 
the auto-correlation of the fluctuating force. Within 
MCT it is assumed Q that its dominant contribution 
arises from pair-density modes p q pk- q - A projection 
of the fluctuating force onto this pair subspace gives 
7fc(*) ~ £,,«,< V* k _ q V ql ,k- q >F q , k - q . q >.k-q>(t) where V p , q 
represents static projections of the fluctuating force and 

Fp,q,p',q'{t) = (P-p(Q)P-q{Q)Pp'(t)Pq'(t)) is a four-point 
correlator with projected dynamics. On an ad hoc ba- 
sis, static and dynamic correlations are then subjected 
to Gaussian factorization. An additional convolution ap- 
proximation for the statics reduces the memory func- 
tion to 7fc(t) ks J dqw q M- q F q (t)F\ k _ q \(t), with a weight- 
factor Wp yq containing only static information. Under 
these approximations for 7fe(i) equation is closed in 
Fk(t) and we have arrived at standard MCT. 



Various versions of extended mode coupling theories 
(eMCT) have been put forward to account for activated 
processes . These theories perturbatively invoke a cou- 
pling to current modes. But recent simulations of Brown- 
ian systems @ , where the momentum current should not 
effectively couple to slow relaxation, have shown that de- 
viations from MCT can be as large as they are in New- 
tonian systems. Also, recent theoretical work has cast 
grave doubt on the role of current modes in eMCT 0. 
It would appear that perturbative coupling to currents 
cannot provide the local physics necessary to restore er- 
godicity deep within the activated regime. 

Within a theory entirely based on density modes one 
is thus lead to reconsider the applicability of Gaussian 
factorization 0,0,0]. Proceeding in the projection op- 
erator approach one shows that the actual evolution of 
four-point correlations is governed by an equation analo- 
gous to Q that, in turn, couples to six-point correlations, 
and so on. The factorization approximation may thus be 
delayed, a perturbative framework we refer to as gen- 
eralized mode coupling theory (gMCT). Recent imple- 
mentations of gMCT 13, [l0( have demonstrated that the 
inclusion of higher order density correlations appears to 
systematically lower the transition temperature (or raise 
the transition volume fraction), suggesting that some as- 
pect of activated behavior is captured. 

In this letter we illustrate the potential of gMCT to 
account for activated processes, meaning - here and 
throughout - the processes discarded in MCT due to 
factorization. This may be achieved through the non- 
perturbative limit where the full multi-point basis of dy- 
namical correlations is included 11]. Our discussion is 
based on a schematic model as its analytic tractability 
provides deeper insights than could be obtained from 
numerical analysis of more realistic systems. To mo- 
tivate its form we briefly return to Eq. Q and its 
higher order generalizations. It is useful to normal- 
ize 4>k(t) = Fk(t)/Sk- For simplicity and following 
0, 0] we focus on diagonal contributions q = q' in the 
memory function 7fe(£). Treating static correlations as 
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in MCT but retaining the dynamic four-point correla- 
tions </>k u h a (t) = F klMMM (t) / F klMMM (Q) leads to 
7fe(*) ~ / d<l A-q,k-q<f>q,k-q(t) where A Pi g = S p S q w p ^ q 
with the same weight-factor w P:q as in MCT. At second- 
lowest order the projection operators yield 0, 

ipkxteit) + Mfei,fe 2 0fei,fe2(*) + (Tfci.fe * 0fci,fc 2 )(*) = 0- 

with Hki,k 2 ~ Mfei + A*fe2- The fluctuating force for pair 
densities is essentially a product of three density modes. 
We perform a corresponding projection and again only 
consider diagonal wave vector terms in r yk 1 ,k 2 (t)- Sim- 
ple treatment of the static projections as before then 

gives Jkite(t) ~ \PktJ dq^qM-q^qM-qMW + Mfc 2 * 

/ dqk qM - q 4> qMM - q {t)]/{ii kl + Hk 2 ) and fa^kufaft) a 
normalized six-point correlation function. 

This procedure may be continued indefinitely and in- 
duces a hierarchy of evolution equations for dynamical 
multi-point correlations. Our schematic model is ob- 
tained by dropping the wave-vector indices. We replace 
4>k{i) <f>i(t), 0fe 1 .fc 2 02 (t), etc., and since we effec- 
tively do not discriminate between different wave- vectors 
Mfc, + ■ • • + fJ.k„ | — * fin. The memory-functions naturally 
[l2| become 7fe 1 ,...,fe„(t) l— * M^n+lC*)- Neglecting iner- 
tial effects and including the bare viscosities arising from 
the short-time portion of the memory functions we arrive 
at the schematic hierarchy 

4> n (t) + fJ,n(p n (t) + fiA(<j) n+1 * 4> n )[t) = 0. (2) 

From now on we set fj, = 1 as it can be absorbed by 
a rescaling of time. Note that Eq. (j2J) reduces to the F 2 
model when closed with ^(t) — 0?(i) El El- The drastic 
simplifications inherent in our schematic hierarchy clearly 
prohibit a quantitative interpretation of its predictions. 
Nevertheless, one would expect that Eq. (J2J captures the 
general phenomenology of gMCT in a similar way as the 
F 2 model does for MCT. 

The entire analysis of our schematic hierarchy is based 
on one central identity. It applies in Laplace transformed 
representation, 0„(s) = C{c/) n (t)}, where Eq. (0) reads 



4>n{s) 



l + A0 n+ i(s) 



(3) 



Recursion of this equation produces a continued fraction. 
Using various transformations (all mathematical details 
will be given in |l4|L however, we derive an identity that 
directly relates any two functions <f> m and <p n : 



fi m (s) = n„(s) with Q„(s) = r(n) 



Vn(s) 
u n (s)'~ 



(4) 



where T(z) denotes the Gamma function and 



u„{s) = *(s + l,s + n + l;A) + 0„(s)*(s,s + n;A), 
v n {s) = $(s + 1, s + n+ 1; A) - 0„(s)$(s, s + n; A). 



The latter expressions contain the (regularized) confluent 
hypergeometric functions $(a, b;z) = i-Fi(a, b; z)/T(b) 
and ^(a,b;z) 0. Note that Q m (s) = Q n (s) for any 
m,rt implies that f2„(s) is independent of n and thus in- 
variant under @. So any single m (s) determines the 
invariant fi(s) = Q m (s) and from that, in turn, all (f> n {s) 
follow by simply rearranging f2(s) = f2„(s). 

In gMCT the hierarchy (0) is closed through factoriza- 
tion at some order N > 2. In general this amounts to 

(j> N (t) = 0™ 1 O)0™ 2 2 (*)••■ C* (*)■ Here orders ™* appear 
with multiplicity nij > 1 such that AT = nim\ + n^mi + 
... + n k m k . In schematic MCT <0 2 = 0? while gMCT clo- 
sures considered previously correspond to 4>s = <pi4>2 
and </>4 = 0i02 10]. It turns out that for tmj/ factoriza- 
tion there is a critical coupling A c above which the <p n (t) 
do not relax fully, i.e., q n — 4> n (t — > oo) > 0. This can 
be deduced from Eq. associated with each closure 
is a polynomial, say in qi, whose real roots are dynami- 
cal fixed points. Some examples of critical couplings A c , 
where real roots appear, and the corresponding plateau 
heights qi for low order closures are listed in Tab. [I] We 
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TABLE I: Critical couplings A c and associated plateaus qi 
(j>i(t —* oo) for all factorizations of 0jv(£) with N — 2, 3, 4. 



observe that within each order N the value of A c is low- 
est for the factorization 4>n = 4>\ while it is largest for 
4>n = 4>n4>N-n with n the integer part of N/2. More im- 
portantly there is an overall increase in A c as we raise N. 
Indeed, for the factorization 0jv = <j>\ it can be shown 
that A c ~ N for N 3> 1. Including higher order correla- 
tions does not merely affect A c quantitatively but in fact 
allows us to increase its value unboundedly. Implications 
of this result will become clear in the following. 

The dynamics of gMCT under higher order factor- 
ization closures, which has never been investigated, can 
be obtained (schematically) by numerical integration of 
Eq. (0). Due to factorization of 0jv(*) the hierarchy re- 
duces to a set of coupled non-linear integro-diffcrcntial 
equations in {<fii(t), . . . , (j>]y—r(t)\. It is efficiently inte- 
grated by the algorithm of [16|. The results in gMCT 
with N > 3 turn out to be remarkably similar to those of 
MCT. For any given factorization, 4>i(t) displays two- 
step relaxation when A approaches the relevant criti- 
cal value A c from below. As A — > A c the a-relaxation 
time t diverges and 0i(f) assumes the non-zero limit 
?i = <j>i[t —* oo). We have performed a scaling analysis, 
which is based on Eq. (J3J and 4>n(s) = C {J^j 0™/(*)}; to 
precisely characterize the critical properties near A c . It 
can be shown that the scaling equations for our hierarchy 
are in fact formally identical to those for usual schematic 
models t l|. Hence /3-relaxation follows the power-laws 
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FIG. 1: [color online] Exact solutions of the infinite hierarchy, 
Eq. 0, at A = 1, 2, . . . , 10 (full curves), MCT with closure 



= 4>\ at A = 1, 2 (dashed) and gMCT under (f> 7 = 
A = 4, 5, 6 (dashed-dotted). See text for discussion. 



at 



4>i(t) ~ Qi ~ t~ a an< i Qi ~ 0i(*) ~ ^ b i n the early and 
late regimes, respectively, and the a-relaxation time di- 
verges like t ~ 1/ (A c — A) 7 . Strikingly the exponents a, 6 
and 7 for these scaling laws are identical for all factoriza- 
tion closures 4>n — Yli ■ Their values exactly match 
those of the model that is a w 0.395, 6=1 and 
7 = j- + 55, regardless of closure. 

While tempting, the above results cannot be inter- 
preted as supporting the standard MCT picture. The 
occurrence of the robust critical scaling under higher or- 
der closures is systematically deferred along with the crit- 
ical point A c ~ N itself. We turn to the non-perturbative 
limit N — > 00 of avoided factorization to develop a deeper 
understanding. According to our general result (@J all 
functions 4> n {s) are known once we have determined the 
invariant fi(s). A priori the latter depends on how we 
close the hierarchy since we may write ft(s) = fijv(s). 
However, analysis of f2/v(s) shows that for all physically 
reasonable closures the invariant of the infinite hierar- 
chy vanishes fl(s) = lirn/v^oo &n(s) = 0. Thus the 
infinite hierarchy has a unique solution. Rearranging 
^n(s) = O(s) =0 then produces the exact result 



4>n{s) 



$(s + l,s- 



1;A) 



<&(s, s + n; A) 



(5) 



The regularized confluent hypergeometric functions $ are 
analytic in s. Therefore the only singularities in Eq. JSJ) 
are (first order) poles at points {si} where the denomina- 
tor &(si, Si+n; A) = vanishes. The relaxation spectrum 
{.Si} is infinite, discrete and contained in the negative 
real axis. Consequently the inverse Laplace transform 
of Eq. (0 is of the form <p n (t) — r i eSit with r* the 

residues at s = s,. Numerical evaluation of both {si} and 
{rj} is straightforward. Plots of the solutions 4>i(t) of the 
infinite hierarchy are shown in Fig. ^ Upon increasing 
A a single component sq of the spectrum approaches the 
origin. Hence the shape of a-relaxation in our schematic 



hierarchy is exponential with r defined by sot = —t/r. 
From 3>(so, So + 1; A) = one shows that r ~ e A /A for 
A » 1. There is no MCT transition in the infinite hi- 
erarchy at any finite A. Instead, the relaxation time r 
essentially grows exponentially in A. Remarkably this 
hallmark of non-perturbative behavior Q emerges from 
our microscopically motivated dynamical approach. 

We now discuss how the above results merge into a 
consistent picture. In the infinite hierarchy, memory ef- 
fects first appear in the relaxation of cf>i (t) but gradually 
spread into the hierarchy upon increasing A. This causes 
the exponential slowdown in the dynamics of <fii(t). The 
MCT closure 02 = 4>\, while appropriate for small A, 
generally overestimates the relaxation time of 4>2(t). As 
shown in Fig. ^ the MCT solution matches that of the 
infinite hierarchy at A = 1 but already deviates at A = 2. 
When A approaches A c = 4 the feedback through the clo- 
sure drives the MCT solution into a spurious transition 
while, in the infinite hierarchy, the relaxation of 4>\(t) 
only develops a shoulder (see Fig. Q] at A = 4) . This is 
also typically observed in experiments and simulations at 
the point where MCT predicts a glass transition 0- If, 
however, the factorization is delayed to higher order cor- 
relations, then so is the appearance of feedback through 
the closure. We refer again to Fig. ^ which also shows 
gMCT solutions under factorization ^7 = 0304. These 
perfectly match the infinite hierarchy even at A = 4, but 
again deviate as A is raised further towards the critical 
value A c « 8.1049 of this closure. Of course the valid- 
ity of gMCT may be extended to any value A by raising 
N . In the non-perturbative limit the transition is de- 
ferred to A c ~ N — > 00 and we arrive at the ergodic 
solution JSJ of the infinite hierarchy. The fact that the 
(j) n (t) derived from Eq. (JSJ do not obey any factorization 
corresponds to underlying non- Gaussian dynamical fluc- 
tuations 0| ■ In microscopic systems these are a natural 
consequence of cooperative, heterogeneous dynamics - 
not fully accounted for in MCT due to Gaussian factor- 
ization. However, microscopic gMCT offers the potential 
to systematically capture non-Gaussian fluctuations and 
thus the underlying processes that cause them. 

Strikingly the solutions of the infinite hierarchy nev- 
ertheless carry some important features consistent with 
MCT. Figure El shows gMCT solutions under factoriza- 
tion at orders 7V = 2,3,...,7at the corresponding criti- 
cal couplings A c . From our discussion of the dynamics in 
gMCT we know that each of these functions approaches 
its relevant plateau q%, see Tab. [H like a power law t~ a 
with the MCT exponent a w 0.395. Now compare these 
critical solutions to those of the infinite hierarchy, Fig. [21 
at the same values of A: apparently there is an excel- 
lent match throughout the early /3-relaxation regime. In 
this sense the hierarchy appears to approach criticality 
with a /3-relaxation exponent a as predicted by MCT. 
Likewise, we could extract the MCT exponent 6=1 
characterizing late /3-relaxation where <p\ (t) drops below 
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FIG. 2: [color online] Main Panel: critical solutions (dashed) 
for closures (f>\, fafa, <p%, 0203, 0§, 0304 at A c ps 4, 4.8, 
5.7, 6.5, 7.3, 8.1, respectively, and solutions of the infinite 
hierarchy at the same A (full lines). Inset: a-relaxation time 
of the infinite hierarchy (dots) for A = 2, 2.5, . . . , 10 and an 
MCT fit (full curve), r flt = A/(A^ - A" 1 ) 7 , with 7 « 1.765 
fixed but A « 0.7 and A c w 10.1 fitted. 



its plateau. A consistent MCT picture would then re- 
quire power-law scaling of the relaxation time r with a 
particular exponent 7. For fitting purposes this power- 
law is usually expressed in T — T c rather than A c — A, 
which is equivalent (asymptotically) close to criticality. 
Noting that 1/A is generally an increasing function of 
T (and thus represents a temperature-like variable) 
we write rgt = ^/(A^ 1 — A^ 1 ) 7 . Most remarkably this 
power law with the MCT exponent 7 w 1.765 indeed 
fits the relaxation time of the infinite hierarchy well over 
3 < A < 9 covering, as in simulations Q, about two to 
three decades in r, see inset Fig. [21 Thus, within this 
limited dynamical range the solutions of the infinite hi- 
erarchy exhibit power-law behaviors characterized by the 
MCT exponents a,b and 7. The fitted A c w 10.1, lying 
well above the MCT value A c = 4, is consistent with cor- 
respondingly lower T c extrapolated in experiments and 
simulations 0- However, contrary to MCT these power- 
laws do not reflect asymptotic scaling in the infinite hi- 
erarchy In fact, there is no divergence of the relaxation 
time r ~ e A /A at any finite A, and in particular not at 
A c « 10.1 as extrapolated from the power-law fit. 

The entire phenomenology discussed above emerges 
purely from the structure of our schematic hierarchy J5J). 
The qualitative differences between the non-perturbative 
and MCT solutions of Eq. J2J m several aspects resemble 
those found when contrasting experimental or simulation 
data |i|6j,ll7J with microscopic (fc-dependent) MCT, e.g., 
predicted versus fitted T c or the apparent power-law be- 
haviors and deviations away from them. Our analysis 
resolves these discrepancies, leading to a coherent pic- 
ture based upon non-Gaussian dynamical fluctuations. 
To describe the relaxation of supercooled liquids in full 



detail, including, e.g., stretched exponentials, numerical 
analysis of microscopic (fc-dependent) gMCT will be nec- 
essary. One would hope that such an approach preserves 
the successful features of MCT, for instance, accurate 
prediction of the plateau values. Beyond that, however, 
it might not only offer a possibility to systematically 
improve the predictions of MCT, particularly through 
the non-perturbative limit, but also provide a kinetic ap- 
proach to capture dynamical heterogeneity. The impli- 
cations of the work presented here on related issues, such 
as the breakdown of the Stokes-Einstein relation, will be 
explored in future work. 
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